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/■^v , We study the mechanism of transverse energy (Ex) production in Au+Au collisions 

at RHIC. The time evolution starting from the initial energy loss to the final Et 
production is closely examined in transport models. The relationship between 

£^ ' the experimentally measured Et distribution and the maximum energy density 

^ ' achieved is discussed. 
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1 Introduction 



The Relativistic Heavy Ion Collider (RHIC) will provide us with a unique 
opportunity to study matter under extreme conditions. The transition from 
hadronic degrees of freedom to quark and gluon degrees of freedom is expected 
from collisions at these energies. The transverse energy spectrum will be 
among the very first results from RHIC. One important question we would 
like to answer from the day-one physics at RHIC is the -maximum energy 
OS ■ density reached in central Au+Au collisions. Lattice QCD El calculations give 

us the critical energy density for the transition to quark-gluon plasma. In this 
talk, we explore the relationship between the transverse energy spectrum and 
the energy density, and discuss the possibility of experimentally determine 

q ■ the maximum energy density. 

^ \ A general framework for computing and analysing energy-momentum ten- 

sor is introduced. As the first part of this study, we construct the energy- 
momentum tensor for simulated events from transport model RQMD 13. We 
also find the transverse energy, the equation of state, and the energy density 
r> \ for these events. Currently we are using our techniques to study the parton 

j^ ■ cascade model VNI a as well as other event generators. 

2 Energy-Momentum Tensor in Transport Models 

In transport model, the energy-momentum tensor 

T^(g,t) = YlJ d 3 p^Mx,p,t) , (1) 
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and the particle number current 



.7' 



'(*>*) ^E/^W'i*'*) 



(2) 



where fi(x,p,i) is the distribution functions for particle type i. 

In a relativistic cascade, each particle is represented by a point in both 
position and momentum space, and the distribution function is an ensemble 
average of the ^-functions, 



Mx,p,t) = (y2S 3 (x-r k (t))S :i (p-p k (t))j . 
The energy-momentum tensor and the particle current become 

T^(x,t) 
and 



I3 k (t)ev „ - 
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(4) 



(5) 



At any given space-time location, T 7 '"', being a symmetric tensor, has 
ten independent elements from which we can obtain ten physical quantities: 
local energy density e, local pressures Vx.Vi.Vz, the flow velocity Vf, and the 
orientation of the principal momentum axises. Following the convention of 
Landau and Lifshitz □, the local rest frame is defined as the frame in which 
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The Lorentz boost of T^ v to the above form gives us the flow velocity. The 
momentum tensor T lJ canJpe diagonalized by performing a rotation. After 
the boost and the rotation I 



rpf_iV 



The local particle density 
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Figure 1. (a) Transverse energy at mid-rapidity as a function of time for central (b = 0) 
Au+Au collisions at RHIC from RQMD v2.4; (b) Local energy density, denned by Eq.(6), 
as a function of time at various positions (z = and r t = 0, 4, 8fm) for the same set of 
collisions. 



where u^ = (7,71?) is the velocity four vector. We further define 
V = {Vi + V 2 + P 3 )/3 and T = V/p . 



(9) 



If there is local thermal equilibrium, then V would be the pressure and T 
would correspond to the temperature. 

3 Event Generator Studies 

There are, in general, two stages in the evolution of E t . Initially, E t increases 
due to the transverse excitations from interactions. Afterwards, the longitudi- 
nal expansion results in a decreasing E t . What we measure in the experiment 
is the final E t at the end of this expansion. The details of the E t evolution, 
like all other observables in the relativistic heavy ion collisions, is model de- 
pendent. In order to get a reasonable estimate of the energy density through 
measurement, we need a systematic study of all plausible models. The trans- 
verse energy E t and the energy density e from the transport model RQMD is 
shown in Figure 1. 
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Figure 2. For the same set of collisions from RQMD v2.4 as in Figure 1 :(a) The ratio 
of local pressure, defined by Eq. (9), to local energy density as a function of time at various 
positions (z = and r t = 0,4, 8fm); (b) Flow velocities for tt,K, and for other particles as 
a function of time at z = and rt = 6fm. 



A total of 144 central (b = 0) Au+Au events, with full event histories, at 
RHIC energy (100 GeV/A+100 GeV/A) was generated using RQMD version 
2.4 0. The energy- momentum tensor and the particle current is computed 
according to Eqs. (4) and (5), where the small volume V is chosen to be a 
sphere of an 1 fm radius. The local energy density e, pressure V, and the flow 
velocity Vf are obtained from Eqs. (6)-(9). Figure 2 shows the variation of 
V/e and the development of transverse flow velocity. 

In the RQMD events, during the early stage, the longitudinal pressure is 
higher than the transverse pressure. After about 7 fm/c, the system is ap- 
proximately isotropic, but it is still not in a local thermal equilibrium. Figure 
2(b) displays one such nonthermal behaviour in the flow velocity distribution. 
Particles with lighter masses, e.g. pions and Kaons, have greater collective 
velocities than heavier particles, i.e. there is no common flow velocity. 

Clearly, the system in Figures 1 and 2, can not be described by a simple 
Bjorken expansion picture. There are expansions along both the longitudinal 
and the transverse directions, and the system is not in a local thermal equilib- 
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rium. The time evolution of e and E t , in Figure 1, are not exactly correlated. 
So the measurement of E t by itself is not sufficient to determine e. 

The initial state in RQMD is not partonic, so it may not be applicable for 
RHIC. But the model may still have a reasonable parametrisation for stopping 
and for initial transverse energy production, and it has a well tested hadronic 
final state after-burner. Models like VNI, HIJING, VENUS, and UrQMD 
have quite different descriptions of the initial stages at RHIC. By coupling 
these initial conditions with a common hadronic after-burner, we will have 
a better understanding of the relationship between the observables and the 
maximum energy density. 

4 Remarks 

An estimate of the maximum energy density might be possible if,in addition 
to transverse energy, we also have a good measurement of several collective 
observables such as the radial and elliptic flows, the system size from HBT, 
and other correlations . This kind of estimate would still be model dependent, 
so we need to make a survey of all models. 
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